A Convenient Alternative for Series Manipulation via the Translation Operator 
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We derive and discuss a technique for manipulating power series which is complementary to 
standard procedures. We begin with the translation operator, but we express the operator as 
an infinite product instead of expanding it as a series and we apply combinatorial arguments to 
generate the terms in the series in an efficient manner with a minimum of clutter and intermediate 
calculations. The method is effective for developing multivariate expansions, and may also be used 
to manipulate series, e.g. in operations where one must take the reciprocal of a power series or raise 
it to a power that may be fractional or irrational. In the case of two component perturbations, 
we obtain analytic expressions for the expansion coefficients. We use our technique to generate an 
electrostatic multipole expansion as a demonstration of its utility in producing coefficients of special 
functions such as the Legendre and Hermite polynomials. 

PACS numbers: PACS numbers: 02.30.Lt, 02.30.Mv, 02.30.Gp 



I. INTRODUCTION 

Taylor series expansions are an important tool in the 
analysis of the effect of small perturbations where the de- 
velopment of a power series in the perturbing influence 
conveniently organizes successive corrections to the ze- 
roth order result. If a single variable function is written 
as /(xo + Ax), where Aa; is the perturbation of small 
displacement, then the Taylor has the form 

f{xo + Ax) = f{xo) + Ax^f{xo) + l^Ax'^f{xo) + ... 

(1) 

A salient example familiar to the modern physics student 
is the case of the Lorentz factor central to formulae from 
special relativity. The Lorentz factor appears in the ex- 
pression for the total energy of a body in motion, jmoc^, 
where 7 = (1 — mo is the rest mass, and c is 

the speed of light; (3 = v/c is the object's speed relative 
to the speed of light. With a single variable Taylor se- 
ries expansion, we readily obtain the energy of the mass 
when it is at rest, the classical kinetic energy term, and 
the flrst relativistic correction. One has 



E = jmoc^ = (1 - /52)-i/2moc2 



(2) 



For ordinary speeds (i.e. for very small the Lorentz 
factor can be expanded as 

Thus, we flnd for the total energy of the moving mass 
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E = ^moc^ — ma(? + -too'I'^ + ^"^o^'^'/c^ (4) 



(5) 



with the kinetic energy being augmented by the factor 
l + |/3^.The latter correction becomes more important at 
relativistic speeds where the kinetic energy scales more 
rapidly than the Newtonian quadratic dependence. 



For situations involving more than one perturbing in- 
fluence such as the function /(a;o + Ax, + Ay, zq + Az) 
where there are several displacements, a multivariable 
Taylor series expansion is often appropriate. As an ex- 
ample, one could have an expression raised to a fractional 
power such as (1 -I- [Ax + Ay + Az])" where the pertur- 
bations are combined in a single term; in fact we will 
see that broadly speaking, the analysis of expressions of 
the form /(I + [Ax + Ay + Az]) will be facilitated with 
the techniques we report on here. To develop the power 
series in the displacements, one could begin with 



/[l + (Aa; + Ay + Az)] 
+ l^/'{l)[Ax- 



= 1 + / (l)[Ax. 
Ay + Azf + . 



Ay + Az] 



(6) 



From Eq. [6l we see that we have flrst the task of gener- 
ating the primary series and then the evaluation of the 
individual terms, which will be polynomials (i.e. binomi- 
als for the case of two variables, trinomials for the three 
variable case, etc.) raised to successively increasing inte- 
ger powers. 

The case of a perturbation which happens to consist of 
a power series expansion in a small parameter can also 
be a computationally intensive scenario if one wishes to 
calculate a perturbative expansion in the small quantity. 
In the discussion which follows, we will use a; as a generic 
label for the latter. As an example, one may consider 
sec(a;) = 1/ cos(a;), where 

sec(a;) = (1 - x^/2l + x'^/4\ - x^Gl + . . .)"! (7) 
= 1 - i-x^/2l + x^/4l - xVG! + ...)+ (8) 
{~-x^/2\ + x^il - x^Gl + . . .)2 + . . . (9) 

Again, even after the initial expansion, the individual 
terms will need to be evaluated and like powers of x com- 
bined. 

We now discuss a complementary and less laborious ap- 
proach for the development of a Taylor series where the 
perturbation has a composite form such as the multivari- 
able case in Eq. [6] or the situation where the perturbing 
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influence is itself a power series in a smaller parameter, x. 
Our technique is based on the translation operator which 
can induce a coordinate shift by generating a Taylor se- 
ries in one or more variable, and the method eliminates 
intermediate calculations while minimizing clutter. 



II. THE SINGLE VARIABLE CASE 

We comment flrst on the single variable case; excellent 
discussions may be found in standard textsi. Consider 
a function f{xo + Ax) where as before there is a shift 
about xq of Ax. For infinitesimal translations 6x, rea- 
sonably behaved functions can be accurately rendered in 
linearized form with /{xq + Sx) — /{xq) + Sxf (xq). In 
terms of the differential operator d/dx one has, equiva- 
lently, 

fixQ + Sx) = (1 + Sxd/dx)f{x)\^^^^ (10) 
It is possible to build up several shifts with 



f{xo + nSx) = (1 4- 5xd / dx)"^ f {x)\x=x 



(11) 



where there are n copies of the differential translation op- 
erator {\-\-5d/dx). This logic suggests that we may more 
accurately take into a finite shift Ax by combining many 
smaller translations to set up the larger shift. Hence, it 
is reasonable to suggest that 



/(xo -f Ax) = lim (1 -f Ax/N) 



N 



N 



(12) 



becomes an increasingly accurate approximation as N 
becomes very large. The identity lim (1 -t- Ax/N)^ = 

N^oo 

e^^ suggests that we associate a finite translation Ax 
with the action of the exponentiated differential operator 
and well known result 



f{Ax) 



lim (1 -t- Ax/N) 



N 



(13) 



In conjunction with this result, we apply combinato- 
rial logic to develop a series in Axd/dx. Although the 
series will turn out to be nothing other than the stan- 
dard single variable Taylor expansion given in Eq [1] a 
very similar approach will be applied to obtain the terms 
of the multivariable expansion where there are more in- 
tricacies to navigate. Returning to the infinite product 
form, (1 -I- Ax/Nd/dx)'^ , we examine how the power se- 
ries will take shape as N becomes very large. 

The zeroth order term in the power series will clearly 
be unity. To obtain the first order term, we must mix 
in a factor of ^ ^ , and there are a total of N ways of 
doing this. For the second order term, there are N ways 
to choose the first factor of and — 1 ways to 

choose the second factor. However, with the order not 
being relevant, there will be redundant terms, and we 
must divide by a permutation factor of 2! to compensate 



for double counting. Proceeding in a similar way, we 
see that the prefactor for the third order term will be 
N{N — 1){N — 2)/3!. Combining these results, we see 
that the entire series will have the form 



Aa;,^^,, , , ^ Ndf ^ ^N{N-l)d^f 



N 



2\m dx^ 



(14) 

Finally, in the large N limit, factors involving N cancel 
and we obtain the anticipated result 



Aa; d 
'Wdi 



N 



(15) 



III. THE MULTIVARIABLE CASE 

For the multivariable case /(t/i -I- Ai, ?;2 + A2, . 
will effect the translation with 



r(Ai,A2,...,A„ 



exp Ai 



Noting that Ai^ 
becomes 



d_ 

dyi 



d 



+ A, 



dyn 



j, we 

(16) 
(17) 



T(Ai,A2,...,A„) ^cxp 



AyV, this 



V 



(18) 



En route to the multivariable Taylor series, one might 

develop a power series in Ay • V, but the result will be a 
multivariable Taylor series requiring the subsequent eval- 
uation of terms with polynomials raised to integer powers 
as discussed earlier. Instead, just as we did for the single 
variable case, we write the multivariable expansion as an 
infinite product in Heu of an infinite series. The result is 




nAi,A2,. 

A ^ A ^ 

Ai h A2 

dyi dy2 



■ , A„) 



'dy„ 



N 



(19) 
, (20) 



where it is understood that the large N limit will be 
taken. 

Using combinatorial logic in conjunction with Eq. [20] 
will permit the immediate enumeration of specific terms. 
As an example, the coefficient for the linear term A^ 
will be Ai^-^J{yi,y2, ■ ■ - yn) with the TV factors can- 
celing. For quadratic terms, there are two cases. When 
i j there are N{N — 1) ways of collecting terms with 
the factors AiA^; we obtain AiAj g^g^/(yi, 2/2, • • ■) in 
the large TV limit. Note that since i and j are distinct, 
there is no need to divide by a redundancy factor. On 
the other hand, when i — j, we need to compensate 
for double counting by dividing by 2!, and the term is 



rA 



2 d 



^/(yi,y2, . . .y„) as iV 
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In gonoral, one obtains a partial derivative d/dyj for 
each copy of , and we must divide by the permutation 
factor nj\ to compensate for multiple counting. As an 
example, the coefficient of the term A1AIA3A4 would 



be 



:prf{yi,y2,y3,y4,---,yn)- 



21 4! dyidy^dys 

A salient and frequently occurring case where we will 
gain an advantage arc situations in which the perturbing 
elements are combined together in a single term, i.e. as 
in /(I + [Ai + A2 + . . . + A„]). An important example 
is the case of an expression raised to a power a which 
needn't be an integer (though there are nontrivial cases 
that are often of interest such as the reciprocal, where 
a = — 1). Our expression is hence 

/(yi + Ai,y2 + A2,...,y„ + A„) (21) 
= (1 + [Ai + A2 + . . . + A„])" (22) 

What will work in our favor and facilitate the calculation 
both in this special case and the broader class of func- 
tions /(l-|-[Ai-|-A2 + . . .+A„]) is the fact that the partial 
derivatives depend only on the total order of the deriva- 
tive, not the variable with which one differentiates. As an 
example, the derivatives QyJy^Qy^ fiyi,y2, ys, ■ ■ ■ , Vn) = 
^/(2/i,2/2,y3, • • • ,2/n) = a(Q; - l)(a - 2) are identical. 
Let us examine the consequences of this simplifying prop- 



erty. 



For linear terms, the partial derivative ^ yields a 
factor a (independent of the value of i), and this value 
is the expansion coefficient. For the quadratic terms, the 
partial differentiation will yield for each {i,j} the same 
factor a{a — 1). In addition, we will need to divide by 
the permutation factor 2! if the two indices i and j are 
identical. 

In general, as we continue to the third and higher order 
terms, the only things to keep in mind are the overall 
order of the term (the value of the derivative of / is fixed 
in this way) , and the instances of repeated indices where 
it will be necessary to divide by a multiplicity factor for 
each case where there is a repeated index. For brevity, 
we represent the expansion series as (1 + [Ai + A2 + . . . -I- 

^n]r = 1 + [AlAi + A2A2 + ... + AnAn] + [AnAj + 

J412A1A2] -|- . . . A few examples are 
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A12 = a{a - 1); An = ^a(a - 1) 
^1222 = ^a{a - l)(a - 2){a - 3) 



(23) 

(24) 
(25) 



Ai234 = a{a-l){a-2){a-3) 

Our results for the multivariable case also have rele- 
vance for perturbations which are themselves power se- 
ries in a single variable. One may have, for example, an 
expression of the form 



/(.t) = (1 + [cix + C2X^ + c^x^ + . . .])'■ 



(26) 



where the series may terminate or instead be an infinite 
power series. To make use of the machinery we have dis- 
cussed in the case of more than one variable, we sub- 
stitute a Aj for each term in the series (e.g. we set 



Ai = CiX, A2 = C2x'^ , etc). For convenience in sort- 
ing the terms, we choose the index i to correspond to 
the degree of the term it represents. As an example, in 
developing the power series for sec(a;), we begin with 

sec{x) = + (-xV2! + .xV4! - x^ + . . .)], (27) 

and for practical purposes that will become apparent we 
follow the pattern mentioned above and use the sub- 
stitution scc(a;) = 1/(1 -f A2 + A4 -t- Ae -|- . . .), where 
A2 = -x^/2\, A4 = x''/4!, and so on. 

With the formalism we have discussed, we may now 
systematically and conveniently write down the terms in 
the expansion for [1 -I- (Ai -I- A2 -I- . . .)]" , though we men- 
tion one caveat; it is important to list all possible combi- 
nations of the Ai factors in developing the expansion. For 
terms up to fourth order, this is easy to do by inspection. 
For the higher order cases, one can systematically enu- 
merate all of the Aj combinations by listing all possible 
partitions of the integer corresponding to the order, and 
a recursive procedure for the partitioning is described in 
the Appendix. Since the value of the partial derivatives 
depends only on the total number of differentiations, we 
will use the notation d-i where, e.g., di = df /dyi and 
d2 = d^f/dy^ = f /dyidyj Following the rules derived 
earlier, the coefficients of powers of x are 



rfiAi 



acix 



for the first order contribution, 



diA2 



1 

2!' 



:doA 



2^1 



ac2 + —^a{a - l)c\ 



for the second order piece. 



rfiAs + dsAaAi + j^d^Al = 
ac3 + a{a - l)c2Ci + -^"(^ - l)(a - 2)cf 



(28) 



(29) 



for the third order terms. 



diA4 + d2A3Ai + ^d2Al+ 



1 

2!' 



(30) 



irfsAaA? + ^d^Af = 

ac^ + a{a — l)c3Ci + -^a{a — l)c2+ (31) 
^a{a - l)(a - 2)c2cl + ia(a - l)(a - 2){a - 3)4 



for the fourth order contribution. 



IV. THE TWO COMPONENT CASE AND A 
MULTIPOLE EXPANSION 

A. Two component perturbations in general 

We find a particularly simple result for the case of a 
two component perturbation, where in general we have 
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/(1 + [Ai + A2]) Since we will eventually replace Ai + A2 
with cix + C2x'^ (i.e. as occurs often in generating func- 
tions for special functions such as the Legendre and Her- 
mite polynomials) , we will again separate the terms of the 
expansion into groups where the sum of the subscripts is 
the same; in this manner we will automatically have an 
appropriately organized power series in x when we make 
the substitution yi = Cix, y2 = c^x^ . As is discussed in 
the Appendix, the fact that in this context one partitions 
integers into I's and 2's with no higher integers appearing 
greatly Hmits the number of terms which will be gener- 
ated for each order, and the number of terms associated 
with each order rises only linearly with the order n. 

With these ideas in mind, we see that the first order 
contribution is diAi. The second order piece, with two 
terms, is (di^2 + ^<^2Af). The third order component 
also has two terms since there is no A3 term in the pertur- 
bation, and we obtain d2AiA2 + ^dsAf ; Hkewise, there 
being no A4 term, the fourth order piece has the form 
dil^ -I- ^(i3A2Af -I- ^(i4Af, and one continues in this 
manner for the higher order terms. If we specialize to 
the case y\ = CiX and 2/2 = C2x'^, we find for the power 
series in x 

/(I + C,X + C2X^) = /(I) + { ^d^c^c\x+ 
+ (Y^n'^24c} + ^d3C^C?) X^} 

+ {-^d2clc\ + ^d.clcj + ^d,4ct) + 

(33) 

where 0!, 1!, 0°, and c" have been inserted strategically 
and odd and even terms separated to emphasize general 
trends. To obtain expressions for general even and odd 
powers, we combine the terms in parentheses into sum- 
mations and thereby make the dependence on the order 
n more transparent. We deduce for the coefficients of x 
the expressions 

n 

n ^ 

^r=Y. ^^_^,^_^^ r.,.-.<^'<^-^ (35) 

where we write for the full expansion /(l-|-[ciX-|-C2a;^]) = 
/(I) + Aix + A2X^ + A^x^ ^ ... 

B. Multipole expansion for an electrostatic 
potential 

A salient example from electrostatics where we may 
apply this logic is the multipole expansion of the electric 
potential of a charge displaced from the origin by Az as 



shown in Figure 1. If we operate in the polar system, 
the combined potential due to the negative and positive 
charges is 

V{r,e) = -g/r[(l -2cos6'A/r + A^r^)"^/^ 

-(l + 2cosM/r + AVr2)-^/^] (36) 

For distances large relative to the size of the charge ar- 
rangement (i.e. where r ^ A), we may expand in the 
small parameter A/r. 

The perturbation is a two term series in A/r; following 
the procedure described in this report, we set (1 + Ai + 
A2)-i/2 s (l-2A/rcos6'-hAVr2)-i/2 for the first term 
in the electrostatic potential. The analysis of the second 
term, except for a sign difference, proceeds in an identical 
way. Setting Ci = — 2cos0, C2 = 1, and "a;"= A/r, we 
see that the first order contribution will be 

dicix = -(l/2)(-2)cos6l(A/r) = cos6l(A/r), (37) 

and this is Pi(cos 0)(A/r) where the Pi{x) are Legendre 
polynomials of order I. The second order piece is 

(^diC2 + -^d2clj x^ = 

[-1/2 + (l/2)(-l/2)(-3/2)4cos2 e]{^/rf 
= (-l/2-h3/2cos2 6')(A/r)2 =P2(cos6')(A/r)2 (38) 

The third order contribution is 

(^2CiC2 + ^cfsc'i^ = 

[(-l/2)(-3/2)(-2cos^?) + 

(l/6)(-l/2)(-3/2)(-5/2)(-8cos3 0)](A/r)3 = 

(-3/2COS6' + 5/2cos^6i)(A/r)^ = P3(cos6')(A/r)^ 

(39) 

The fourth and higher order terms are obtained in a sim- 
ilar manner. We have for the combined potential 

V{r, e) = COS e + (-3 cos 61 5 cos^ e) + ... 

(40) 

We recognize the first nonzero term as the dipolar con- 
tribution where the factor 2A(5 is precisely the dipole 
moment appropriate to the geometry and the equal and 
opposite charges +Q and —Q. 

C. Additional Examples: Hermite Polynomials and 
the Bernoulli numbers 

Generating functions often appear in the form /[I + 
{cix -t- C2x'^y\, and we may therefore apply the formalism 
used in the context of the electric multipole expansion 
to generate the first several polynomials in such cases. 
For the Hermite polynomials, the generating function 
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gta: t /2 _ gf^^^ -g particularly amenable to our expan- 
sion technique since the derivatives di = d2 = ■ ■ ■ = dn = 
1 are each equal to unity. With e*"^"* — X^i^o 
with Hn{x) = n!a„, we find that Hq{x) ~ 1 and Hi{x) — 
(l!)(ci) = X, while H2[x) and H^{x) are given by 

(41) 

and 

i?3(a:) = 3! [c2Ci + ic?^ = 6 (^-| + = x3 - Sx. 

(42) 



V. APPENDIX: THE PARTITIONING OF 
INTEGERS 

The task of exhausting all possible combinations 
of Aj/i becomes increasingly subtle with increas- 
ing order. As an example, for the second order 
case one has only {A2,A^}, whereas for fifth or- 
der the set of possibilities is much larger, namely 
{As, A4A1, A3A2, AgAf, AiAi, A2A3, Af }. Neverthe- 
less, one can easily list all possibilities for small to moder- 
ate order by following a recursive partitioning procedure 
which we discuss here. We first examine the general case 
where the perturbing influence may have a large (or even 
infinite) number of terms. We then find a much simpler 
result for the case of a two component perturbation. 



Our technique can be used to conveniently obtain the 
Bernoulli numbers Bn from^ 



2^2 



It is useful to write 



(43) 



(44) 



{l+x/2\ + x'^/i\ + ...) 
We see that Bi — dici ~ — i. For B2, we have 



B2 = 2\\ diC2 + ^d2cn = 2(-l/6 + 1/4) = i (45) 



We find for B-i 



B3 = 3! {^diC3 + d2CiC2 + ^dscl 

= 6(-l/24+ 1/6+ 1/8) = (46) 



D. Conclusions 

In conclusion, we have discussed a technique for facil- 
itating the manipulation of power series which relies on 
the structure of the multivariable translation operator 
when written as an infinite product. The technique elim- 
inates intermediate calculations while minimizing clutter. 
The approach has been demonstrated for series raised to 
powers, but the technique may be applied in more gen- 
eral situations where one seeks to expand a function of 
the form /[I + {aix + 02^^ + ...)] where the perturb- 
ing infiuence is itself a power series in a small quantity 
X. In the case of a two component perturbation (e.g. 
as in generating functions for special functions such as 
the Legendre and Hermite polynomials), the process of 
generating the power series is relatively simple, and our 
technique facilitates the derivation of a analytical closed 
form expression for the coeflicients of each order of x. 



A. The general case 

Let us denote the set of all possible partitions of the 
integer n with curly braces, (i.e. as {n}). It is easy to 
see that for n — 1, there is only {1} — [I]. 

Next, we discuss a systematic procedure for generat- 
ing partitions by examining successively higher integers 
beginning with n — 2 and continuing through n — 6. 
One always begins with the unary partitioning [2], just 
the number itself enclosed in square brackets. The next 
step is to generate all unique binary partitions, and for 
n = 2 there is one of these, [1, 1]. To represent the set of 
possible partitions, we use the group addition symbol © 
to indicate the combination of the individual partitions. 
E.G., for n = 2, we have 



{2} = [2] ©[1,1] 



(47) 



For n = 3, we proceed in a similar way, again, form- 
ing all possible binary partitions and finding {3} — 
[3] © [{2} , 1]. Now, substituting the partitioning {2} de- 
termined for n — 2, we obtain 

{3} = [3]©[[2]©[1,1],1] (48) 

Exploiting linearity and distributing yields 

{3} = [3] ©[2,1] ©[1,1,1] (49) 

for the partitioning corresponding to n = 3. 
For n = 4, we have 

{4} = [4] ©[{3},!]© [{2}, {2}] (50) 

Acting recursively and inserting the relations we obtained 
for {3} and {2} yields 



{4} = [4]©[[3]©[2,1]©[1,1,1],1] 
©[[2] ©[1,1], [2] ©[1,1]] 



(51) 



Distributing and eliminating redundant terms yields 

{4} = [4] © [3, 1] © [2, 2] © [2, 1, 1] © [1, 1, 1, 1] (52) 
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We may proceed in the same manner for n = 5 and 
higher, but the structure of the partitionings {2}, {3}, 
and {4} suggests a more efficient and straightforward 
approach akin to counting. To emphasize the counting 
character of the partitioning process, we have taken care 
to rank the specific partitions in decreasing "dictionary" 
order, where the first entry in the partitioning determines 
priority and ties are broken by subsequent entries where 
needed. We illustrate the systematic counting process for 
the case n = 6, and we obtain {6} in this manner with 
no discarded terms. 

The highest possible number which may appear in {6} 
is just a single "6", so the first entry in the partition- 
ing is [6]. The next highest entry would be a "5", and 
the only possible partition is [5,1]. The highest rank- 
ing entry in the next tier is [4, 2] , and we may decre- 
ment by breaking up the lowest reducible entry; we 
hence split the "2" into two "l"'s yielding [4,1,1]. For 
the next level, the highest ranking entry is [3,3]. The 
only way to generate a lower entry is to reduce the sec- 
ond "3", yielding [3,2,1], followed by [3,1,1,1]. The 
next level and the final tier will consist of the entries 
[2, 2, 2] ffi [2, 2. 1, 1] e [2. 1, 1. 1, 1] © [1, 1, 1, 1, 1, 1]. To fur- 
ther illustrate the decrementation method, we generate 
the partitioning for n = 8 by counting down from [8]. 

{8} = [8] ©[7,1] 
©[6, 2] © [6, 1, 1] 

©[5, 3] ©[5, 2,1] ©[5, 1,1,1] 
©[4, 4] © [4, 3, 1] © [4, 2, 2] © [4, 2, 1, 1] © [4, 1, 1, 1, 1] 
©[3, 3, 2] © [3, 2, 2, 1] © [3, 2, 1, 1, 1] © [3, 1, 1, 1, 1, 1] 
®[2, 2, 2, 2] © [2, 2, 2, 1, 1] © [2, 2, 1, 1, 1, 1] © [2, 1, 1, 1, 1, 1, 1] 
©[1,1,1,1,1,1,1,1] (53) 



B. The case of a two component perturbation 



With a two component perturbation such as A1-I-A2 = 
Cix + C2X^, the process of partitioning the integers for 
each order in x simplifies considerably, and can be done 
essentially by inspection for any order. As an example, 
we have for the first several orders 



{3} 



{1} = [1]; {2} =[2] ©[1,1] 
[2, 1] © [1, 1, 1]; {4} = [2, 2] © [2, 1, 1] © [1, 1, 1, 1] 
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Figure 1: The charge arrangement with Q and —Q on the z 
axis and separated by a distance 2A. 



